Fluctuations and scaling in creep deformation 
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The spatial fluctuations of deformation are studied in creep in the Andrade's power-law and the logarithmic 
phases, using paper samples. Measurements by the Digital Image Correlation technique show that the relative 
strength of the strain rate fluctuations increases with time, in both creep regimes. In the Andrade creep phase 
characterized by a power law decay of the strain rate et ~ t~~ , with 9 fa 0.7, the fluctuations obey Aet ~ £~ 7 , 
with 7 « 0.5. The local deformation follows a data collapse appropriate for an absorbing state/depinning 
transition. Similar behavior is found in a crystal plasticity model, with a jamming or yielding phase transition. 
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The creep of a sample under a constant load exhibits of- 
ten three rheological phases in materials from metals and al- 
loys with a crystalline structure to rocks and composites: the 
initial power-law decrease of the strain rate or primary creep, 
quasi-stationary secondary creep, and finally tertiary creep ap- 
proaching final failure. It is important to understand such em- 
pirical laws of scaling in material deformation. The Andrade 
creep law, stating that the deformation rate decays in time as 
-^ = et ~ t~ e , with 9 « 2/3, originates already from 
1910 ljj]. The phenomenological creep laws |2J] are valid for 
many materials from metals to polymeric materials to com- 
posites like ordinary paper [3] to rocks [4]. In the case of met- 
als the dynamics is related to dislocations and the structures 
they form. It is intriguing why there should be such universal- 
ity, which seems to be "coarse-grained", not dependent on the 
microscopic details. 

In this vein, collective phenomena as fluctuations and 
avalanches Q5[] have been recently studied in crystal plastic- 
ity Bo- lldl . viscoelastic fluids and networks Illl4l3ll . granu- 
lar materials and amorphous glasses [ 14] and in fracture 1 15- 
llal . The dynamics of simple dislocation assemblies exhibits 
a phase transition at a critical external shear stress value <t c , 
similar to the jamming one seen in granular media and glasses 
lil9H2lll . with the order parameter (the deformation rate) ex- 
hibiting a power-law decay in time close to the transition, as 
the Andrade law BH . 

To explore the physics of creep and to see if it could be 
related to such collective behavior, we measure local creep 
rates in primary creep. The main question is: what kind of 
fluctuations does the creep show, and how do they scale with 
time? In non-crystalline materials strain variations implies a 
localization of deformation which is unrelated to dislocation 
activity. So, another issue is: what kind of universality ex- 
ists in creep deformation and in local creep strain rates? By 
digital image correlation analysis, we show that the spatial 
fluctuations of creep strain rates follow a scaling law during 
the Andrade creep. The ratio of fluctuations to mean rate 
increases with time. This feature persists in the logarithmic 
creep regime, where e(t) ~ log(t), or et ~ t -1 . The critical 
properties of the power-law creep are tested by applying the 
theory of non-equilibrium phase transitions. We find that spa- 



tial creep deformation statistics follows its own scaling and a 
scaling function. Similar behavior is encountered in a simple 
discrete dislocation or crystal plasticity model, which exhibits 
Andrade creep for applied stress values close to the critical 
stress a c of the dislocation jamming transition. This stress 
separates an active phase with a constant average creep rate 
from a low-stress phase where the deformation rate dies out 
after a transient [22]. The temperature enters indirectly only 
via the dislocation mobility. The most natural explanation is 
that Andrade creep derives from a non-equilibrium phase tran- 
sition between "jammed", or frozen, and "flowing", or active, 
states 1 19]. Usual materials science constitutive equations and 
creep laws are based on phenomenological equations, usually 
related to an thermally activated process [24]. They may fail 
to cope with correlated, intermittent, dynamics. 

To induce creep in experiments, we applied a constant stress 
to a sheet of paper in a tensile-test geometry (Fig. la). Tests 
were made in humidity and temperature controlled conditions, 
33 ± 1 % of relative humidity and 33 ± 1 °C of tempera- 
ture, with the typical applied loads around 40 N, correspond- 
ing to a stress of 13 MPa. We used ordinary office paper, with 
30 x 100 mm 2 samples loaded in the cross-machine direc- 
tion. The sample is imaged with SensiCam 370KL0562 digi- 
tal camera with a low thermal noise ratio, at 12 bit gray scale 
resolution. The exposure time in measurements was 200 ms, 
and the sampling frequency of the digital imaging was 0.1 Hz. 

In sample scale images, a printed speckle pattern (Fig. lb) 
was used due to its sharp contrast and translational invariance. 
A set of experiments with larger magnification and smaller 
load was also made with an image size of 3x4 mm and using 
specially prepared sheets, with 5 % of the fibres colored for 
contrast. The larger magnification and a smaller load make 
it possible to follow the entire experiment with Digital Image 
Correlation (DIC), beyond the primary creep regime. 

Experimentally, a time-series of displacement fields is de- 
rived using the DIC-technique H25I427I1 (Fig. la). The algo- 
rithm describes both the image and the deformation using the 
B-spline model. The algorithm finds the deformation function 
by using the multiresolution approach for the minimization, so 
that the image and the deformation model is refined every time 
the convergence is reached. The criteria for the convergence 



is the sum of squared differences (SSD), i.e. the deformation 
function is applied to the original image and the SSD is com- 
puted against the deformed image. The algorithm is described 
in Ref. [ 28] . A discrete picture, in space and time, ensues of 
the local creep strain rates |29]. An example (Fig. lb) shows 
the variations of the rates for one At. For the spatial fluctua- 
tions of the local strain rate e\ J , we define an evenly spaced 
grid on an image, which consists of 49x78 points. The stan- 
dard deviation of local strain rates Ae t is computed over the 
grid values on the lower half on the sample where absolute dis- 
placements are large. This is further illustrated in Figure la, 
which shows how the creep fluctuations are measured from 
the experiment. 

A feature visible in experiments is the disordered structure 
of paper, two-dimensional over scales larger than the out-of- 
plane thickness [30]. One should note that the non-uniformity 
results in a locally varying stress field, even in the elastic de- 
formation regime, in analogy to the internal stresses of mate- 
rials with dislocations. 

The Andrade creep results are confronted with those of the 
simulations of a two-dimensional discrete dislocation dynam- 
ics (DDD) model flUKH. The N point-like dislocations 
with Burgers vectors b n = s n bu x (with s n — ±1 and n = 
1 . . . N) parallel to the glide direction can be thought to rep- 
resent a cross-section of a single crystal with straight parallel 
edge dislocations and a single slip geometry [6]. The disloca- 
tions interact through their anisotropic long-range shear stress 
fields, and perform over-damped glide motion only driven by 
an applied external shear stress a ||6|,|23|,|29j]. Fluctuations are 
studied in the vein illustrated by Figure 1 a by measuring the 
local strain rates 
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over a time scale At within sub-volumes of linear size I, with 
x n (t) the position of the nth dislocation within its glide plane. 
Figure 2 depicts the Andrade creep results for the standard 
deviation of the creep rate probability distributions. For both 
paper samples and DDD simulations (performed at a = a c « 
0.025, in dimensionless units 129J) we observe a typical An- 
drade law of the form e t oc t~° J or an Andrade-exponent 9 
close to 2/3. The standard deviation follows also a power-law, 
and scales as 



Ae t oc t 



-0.5 
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so an exponent 7 ~ 0.5 is found. The "second" Andrade 
or fluctuation exponent 7 is slightly higher for experiments 
(0.55), while in the DDD simulations a value of 0.5 is in- 
dicated. Evidently this exponent is smaller than that of the 
mean creep rate. Thus fluctuations increase in relative terms, 
and creep becomes more intermittent with time. The accuracy 
of the DIC method poses a limit to following the power-law 
of Eq. (1), apart from the maximum length of the power-law 
creep regime obtainable by tuning the parameters. In the DDD 



model the scaling range is controlled by the proximity of the 
jamming transition. 

At lower stresses, logarithmic creep is found with no sig- 
nature of the Andrade phase. Figure 3a shows that a fluctua- 
tion scaling can be extracted in this regime. Still, the fluctua- 
tion magnitude decays slower than the creep rate. Their ratio 
might be related in each sample to the transition to tertiary 
creep, which might indicate a zone-like localization of subse- 
quent deformation. The simultaneous measurement of Acous- 
tic Emission indicates no substantial damage accumulation or 
micro-cracking until close to the transition to tertiary creep. 
Thus DIC measures the statistics of viscoplastic deformation, 
and not deformations related to eventual brittle fracture. Mi- 
croscopically in paper, the fibers and the network structure can 
both deform irreversibly. For the sample-to-sample variations, 
we find here a Monkman-Grant relation 11151 13211 between the 
time at which a minimum creep rate occurs in a sample t m 
and the final failure time t c , t m oc 0.83 t c . The relevance of 
fluctuations for the Monkman-Grant law would be interesting 
to understand. 

The theoretical interpretation of the scaling of the creep 
fluctuations can be analyzed within and with the aid of the 
DDD model. One question is how this "second Andrade- 
scaling" is related to the usual primary creep exponent? As- 
suming intermittent, localized creep activity the active frac- 
tion of the regions of the specimen account for the creep 
strain measured in a time step. This argument leads to a self- 
consistency condition between the "number" of such active 
regions as a function of time, and the typical time-dependent 
creep strain rate in each such region J29I1 . If the statistics 
would originate from the decay of the fraction of localized, 
active regions, themselves deforming at a rate independent of 
time, it would indicate a 7 = 1/3, half of the Andrade 6, in- 
stead of the measured 0.5 . . . 0.55. This line of reasoning 
is related to classical, dislocation-based theory arguments of 
Andrade creep [33, 34]. 

The nature of the creep dynamics seems to exhibit univer- 
sality typical of phenomena arising from phase transitions. If 
a wide variety of materials follows such behavior, it calls for 
a continuum description, independent of microscopic details. 
Statistical physics offers concepts such as depinning transi- 
tions of elastic manifolds in random media. These have been 
applied to coarse-grained dislocation plasticity 13511 . and are 
examples of absorbing-state phase transitions. Similarly to 
what is seen in the DDD model, such theories present or pre- 
dict a power-law relaxation of the order parameter close to 
the critical point. The jamming/absorbing state transition sce- 
nario can be further tested by the ansatz that the creep rate is 
the "order parameter" of the phase transition. One collapses 
the total deformation histograms, P(e, t) (scaling out the av- 
erage deformation at time t) as one should do 13611 for the in- 
tegrated order parameter. Both the DIC (Figure 3b) and DDD 
data ([29]) exhibit a reasonable data collapse, indicating crit- 
ical behaviour, close to such a phase transition. Further evi- 
dence for the origin of the fluctuations in collective phenom- 
ena can be sought in the correlation functions C(x — x',t — t') 



of strain fields, calculated as for an "interface field" e(x,t). 
Preliminary results imply, that the temporal scaling yields the 
Andrade exponent II29I1 . The fact that the jamming transition 
picture works for the experiments seems as such surprising 
for a non-crystalline material, though a possible universality 
of the Andrade creep would argue for the opposite. 

It remains to extend the measurements to study the small 
scale spatiotemporal features, such as localized events 13711 . 
The creep fluctuations could be investigated in other sys- 
tems such as two- and three-dimensional colloidal crystals or 
glasses 11381 13911 . Future research directions include consid- 
ering the fluctuations for various load histories, an example 
being fatigue, and the role of "aging" and sample-relaxation. 

Local creep rates exhibit a fluctuation scaling, both in a 
disordered, viscoplastic material and in a model of dislo- 
cation plasticity. This presents a challenge: what kind of 
coarse-grained theory for the deformation field would repro- 
duce those signatures? Such fluctuations are beyond classi- 
cal, phenomenological creep models. The phase transition 
scenario of Andrade creep relegates the temperature to a sec- 
ondary role, in the DDD model via the viscous dynamics of 
the dislocations. In depinning, for stresses below the critical 
one, thermal fluctuations induce a transient towards a ther- 
mally activated steady-state creep [40], which could relate to 
the logarithmic behavior in experiments. 

Non-equilibrium scenarios include various models for the 
plasticity of amorphous materials fll3l ll4J kll 14211 where the 
presence of localized events has been recognized e.g. in the 
so-called Shear Transformation Zone theory. To summarize, 
the main issue is that creep deformation exhibits spatial fluc- 
tuations which scale in time, indicating that the different phe- 
nomenological creep regimes could be understood in terms of 
phase transitions. 
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FIG. 1. The experimental scenario. From a pair of digital images at 
a time-interval At the local deformations are extracted at a grid of 
N x M points, with N = 49 and M = 39 on the lower half plane 
of the sample. The strain fluctuations are measured via the time- 
dependent standard deviation, and compared to the mean creep rate, 
here in primary/ Andrade creep. In the DDD simulations, the system 
is similarly divided into boxes of linear size I, with local strain rates 
e l t ' J . b) A digital image of a paper sample on a scale of 40 mm 2 . Su- 
perposed is a typical deformation grid for a time difference At of 10 
seconds. The color scale indicates the degree of local creep deforma- 
tion (blue: small, red: large). In the background: the experimental 
setup. The visible speckle pattern is printed and designed to have a 
structure and contrast appropriate for the DIC method. 
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FIG. 2. (Color online) Average creep rate and the standard deviation 
of the local creep rates, a) Typical experimental data. The sample av- 
erage lifetime is 800-1600 seconds. The primary-to-secondary creep 
transition takes place continuously at 150-200 seconds. 16 experi- 
ments are used. The effective Andrade exponent is about 0.7 in the 
range 10 to 50 seconds, and becomes close to one (signalling sec- 
ondary creep) above 150 seconds, b) The Andrade's law and the 
fluctuation scaling for the DDD model for three different box sizes 
to compute the local rates, at a ~ a c . The early-time cross-over in 
the fluctuations is due to computation time-interval of At = 100. 



FIG. 3. (Color online) a) Experimental results from logarithmic 
creep using a smaller imaging area, 5 second image intervals, and 
a smaller load, showing the average creep rate and the fluctuations. 
The former decays faster, and at the end these two become compa- 
rable, on the observation scale. The data is an average over 9 ex- 
periments, b) Data collapse of the primary creep strain distribution 
P(e, t) from the experiments. The local deformation distributions 
can be collapsed using the creep strain as the scaling parameter. 
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